On the stability of the generalized QCD ghost dark energy 
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We investigate the stability of the generalized ghost dark energy model against perturbations in the 
FRW background. For this purpose, we use the squared sound speed v 2 whose sign determines the 
stability of the model. We explore the stability of this model in the presence/absence of interaction 
between dark energy and dark matter in both flat and non-flat geometry. In all cases we find almost 
a same result. Due to the existence of a free parameter in the model, the model is theoretically 
capable to lead a dark energy dominated stable universe. However, observational constraints rule 
out such a chance. In conclusion, we find that the generalized ghost dark energy can not lead to a 
stable universe favored by observations at the present time. 



I. INTRODUCTION 



Seeking a convincing explanation for the present acceleration of the cosmic expansion p|, a variety of models 
have been proposed in the literature. These models should be consistent with observational evidences supporting 
the acceleration phase of the universe expansion 2-5J. All of the models can be categorized in two different groups 
the dark energy (DE) 0-[H and the modified gravity [TH3 models. Both of these approaches are extensively 
investigated in the literature during the past decade. In this paper we are interested in DE approach. DE models 
S^p. 1 are based on the assumption that the correct theory of gravitation in large scale is the Einstein's general relativity 
(GR). In this approach it is assumed that there exist an exotic type of energy which its equation of state parameter 
is negative and push the universe to accelerate. 
• i-h 1 One main problem in solving the DE problem is introducing new degrees of freedom which may lead to inconsisten- 
cies. Thus avoiding such problem we are more interested in DE models, based on already presented degrees of freedom 
in physics. One interesting model of this category is the so called ghost dark energy (GDE) model proposed recently 
|20l - t22| . Seeking a solution to U(l) problem, the so-called Veneziano ghost has been proposed in the low energy 
'— effective QCD where they are completely decoupled from the physical sector [23T - I261 ] . However, they contribute to the 
vacuum energy in curved space or time-dependent background. In [2lj the authors discussed that this vacuum energy 
can play the role of DE. This contribution to the vacuum energy in curved space-time is proportional to Aq CD H, 
\J3 \ where H is the Hubble parameter and Aq C£) is QCD mass scale. With Aqcd ~ lOOAfeF and H ~ 10 _33 eV , Aq CD H 
gives the right order of magnitude ~ (3 x 10~ 3 eV^) 4 for the observed DE density [2l[. Following this interesting 
coincidence, various aspects of GDE have been studied [27rl35l l4fH . 

In [37j |. the authors discussed that the contribution of the Veneziano QCD ghost field to the vacuum energy is 
not exactly of order H and a subleading term H 2 appears due to the fact that the vacuum expectation value of the 
energy-momentum tensor is conserved in isolation |36[. They argued that the vacuum energy of the ghost field can 
be written as H + 0(H 2 ) : where the subleading term H 2 in the GDE model might play a crucial role in the early 
evolution of the universe, acting as the early DE. Based on this idea people considered the role of this version of 
the Veneziano ghost field energy to the DE problem and tried to see if there exist a better agreement between this 
latter model and the observations. We call this model as generalized ghost dark energy (GGDE). In this model the 
energy density is written in the form prj — aH + j3H 2 , where (3 is a constant HU- It was shown [38[ that taking the 
subleading term H 2 into account can give bette r ag reement with observational data compared to the usual GDE. The 
GGDE has attracted a lot interest recently [39l l4dl| It was shown 38, 39], that the subleading term H 2 in the energy 
density has negative contribution compared to the leading term H. This result first pointed out in [38| for fiat FRW 
universe and then confirmed for non-flat universe j39|. 

Based on all mentioned above, in this paper we would like to investigate the effects of the subleading term on the 
background perturbations and stability of the GGDE model. It is worth mentioning that in [28] . we found that the 
GDE is instable against perturbations and a stable DE dominated universe cannot be achieved in such a model. 

From observations we know that our universe is in a DE dominated phase. Thus any viable DE model should result 
a stable DE dominated universe. One simple way to check such a stability for any new DE model is to discuss the 
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behavior of the square sound speed (v 2 = dp/ dp) [41| in a DE dominated universe. The sign of v 2 s plays a crucial role 
in determining the stability of the background evolution. If v 2 < 0, it means that we have the classical instability of 
a given perturbation. In contrast v 2 > 0, leaves chance for greeting a stable universe against perturbations. However, 
this does not enough insight to say the model is surely stable but at least can show sounds of instability in the model. 
This approach has been used for exploring some DE models. For example in [12, IH| the authors investigated the 
behavior of the square sound speed for holographic DE as well as the agegraphic DE models and found both of these 
models are instable against background perturbations. Also it was shown that chaplygin gas and tachyon DE have 
positive squared speeds of sound with, v 2 = —w, and thus they are supposed to be stable against small perturbations 
[44l Hfl . The stability of the original GDE was studied in [28[ , and we found the GDE is not capable to result a 
stable DE dominated universe. The issue was investigated for flat, non-flat in the presence/ absence of the interaction 
between DE and dark matter. Also, a same procedure was considered in [46[ to show the stability of the GDE in the 
chameleon Brans-Dicke theory. 

In this paper our aim is to discuss the chance of the GGDE model in resulting to DE dominated universe. In recent 
years several signals have been detected, implying a small interaction between DE and DM is possible. As an instance, 
observational evidences provided by the galaxy cluster Abell A586 supports the interaction between DE and DM [47| . 
Also evidences from CMB and also supernova measurements of the cubic correction to the luminosity distance favor 
a positively curved universe (4|| [49j|. Due to the above observational evidences we consider the stability of GGDE 
model in the presence of interaction in both flat and non-flat background. 

This paper is organized as follows. In the next section, we review the GGDE model in flat universe and discuss 
its instability against perturbation. In section IIII1 we explore the stability of the GGDE model in the presence of 
the interaction term in flat universe. Section HVl is devoted to stability of the GGDE model in a interacting, non-flat 
universe. We summarize our results in section fVl 



II. STABILITY OF NON-INTERACTING GGDE MODEL 
A. Review of the non-interacting GGDE in flat universe 



At first we briefly discuss the non-interacting GGDE in a flat FRW background. For flat universe, the first 
Fricdmann equation read 

H 2 = ^( Pm + P n), (1) 

where p rn and po are, respectively, the energy densities of pressureless matter and DE. The dimensionless energy 
density parameters arc defined as 

_ 8irGp m _ 8irGp D 

According to these definitions, the first Friedmann equation (jTJ) can be rewritten as 

n m + n D = i. (3) 

The conservation equations also read 

p m + 3Hp m = 0, (4) 

p D + 3H PD (l + w D ) = 0, (5) 

which imply that matter and DE components are separately conserved. 
The energy density of GGDE is defined as [38| 

PD =aH + (3H 2 , (6) 

where a is a constant of order Aq CD and Aqcd is QCD mass scale. In the original GDE with Aqcd ~ lOOMeV and 
H ~ 10~ 33 eV , AQ CD iJ gives the right order of magnitude ~ (3 x 10~ 3 eV) 4 for the observed DE density In the 
GGDE, p is a free parameter which can be adjusted for better agreement with observations. 
Taking the time derivative of relation ((5]), we obtain 

p D = H{a + 2(3H). (7) 
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Also differentiating (Q} with respect to time lead 

H = -4irG PD (l + u + w D ), (8) 
where u = p m / Pd- Replacing this relation in ([7]) and also using the continuity equation ([5]) we get 

(1 + w D )[3H - 47rG(a + 2$H)\ = 4nG(a + 2/3H). (9) 
Solving the above equation for wd and noting that u — jf£ as well as 

— (a + 20H) = -2- + —^, (10) 

we obtain 

where £ — 87T ^P _ j g c i ear that this relation reduce to its respective one in the GDE when £ = j27| . 

It is easy to see that at late times where SId — > 1, — » —1, which implies that the GGDE in a flat universe 
mimics a cosmological constant behavior. Also we find that wd of the GDE model cannot cross the phantom divide 
and the universe has a de Sitter phase at late time. 

Next we turn to the deceleration parameter which is defined as 

ad H . , 

where a is the scale factor. Using Eq. (JSJ) and the definition of £7 £> in © we get 

H 3 

-=2 =--n D (l + u + u;zi). ( 13 ) 

Replacing this relation into (fT2")l . and using (TTTI) we find 

= 1 3 {-Tin ( . 

q 2 2 (fir, +£-2)" [ ' 

One can easily check that the deceleration parameter in GDE is retrieved for £ — |27j. At late time where the DE 
dominates (£Id — > 1), independent of the value of the £, we have g = — 1. Besides, taking ftno — 0.72 and adjusting 
£ = 0.01, obtain qo ~ —0.34 for the present value of the deceleration parameter which is in agreement with recent 
observational data |50f | . 

The evolution equation of the GGDE can be obtained as [3j| 

cMp (1 - ggKg - »£,) , 

dlna 2-fi -£ ' 1 J 



B. Stability of the model according to the square sound speed parameter 

The parameter we use through this paper as a factor of stability of the model is the square sound speed v^. In 
classical theory of perturbation we assume a small fluctuation in the background energy density and we would like to 
see if the perturbation grows or will collapse. In the linear perturbation regime, the perturbed energy density of the 
background can be written as 

p(t,x)=p(t) + Sp(t,x), (16) 

where p(t) is unperturbed background energy density. The energy conservation equation (V^T^" = 0) yields [4l| 

Sp = v 2 V 2 Sp(t,x), (17) 

where v 2 = ^ is the squared of the sound speed. Solutions of equation (|17p include two cases of interest. First when 
v 2 is positive Eq. (fTTJ) becomes an ordinary wave equation whose solutions would be oscillatory waves of the form 
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FIG. 1: Evolution of squared sound speed against fin for noninteracting GDE (dashed curve) and GGDE (continues curve) 
in fiat background for £ = O.f. 



Sp = Spoe iult +' lk - x which indicates a propagation mode for the density perturbations. The second case of interest 
occurs when is negative. In this case the frequency of the oscillations becomes pure imaginary and the density 

perturbations will grow with time as 8p = 8poe u ' tJr x . Thus the growing perturbation with time indicates a possible 
emergency of instabilities in the background. 

Here we would like to obtain the sound speed in a flat FRW background filled with pressureless matter and GGDE 
while matter and DE components are separately conserved. The definition of the sound speed reads [22| 

v 2 s = -j- = - = -w + w, (18) 
dp p p 

where in the last step we have used P = wp. Using the conservation equation ([SJ, we obtain 

(19) 



p 3H(l + w D )' 
Also taking the time derivative of Eq. (fTTj) yields 

dw ■ _ {l D (2£-Sl D )+Z(2-£) i 



%hD = d^ n ° = ng,(2-n g + fl a Qd - (20) 



Replacing Eqs. (fT9l) and (|20|) into (fl~8|) and also using (fT5|) we get 



2_ o (€-nu)(nu-i) 



v '~ 2 (n D + 0(2-^ + 0^ (21) 

where we have used 4z = H d p - . One can easily see that setting £ = 0, this result reduces to its respective relation in 
[28j . Having the vl at hand we are ready to discuss about the stability of perturbations. One can easily see from (p?Tj) 
that Vg is negative provided £ > flo- One should note here that can be positive up to present time if £ is limited 
to suitable values and thus £ plays a crucial roles in the stability of perturbations background. However, according to 
[Hj], the presented range of parameter j3 does not let us to set £ to suitable values. This result indicates that due to 
the negativity of the squared sound speed every small perturbation can grow with time which leads to an instability 
in the universe. Thus we can not expect a noninteracting GGDE dominated universe in the future as the fate of the 
universe. The evolution of versus flo is shown in Fig. m 



III. STABILITY OF INTERACTING GGDE IN FLAT UNIVERSE 



Based on a traditional and historical manner usually people consider dark matter and DE separately. However, 
these days observations detect signals of interaction between DM and DE (42|. Also from theoretical point of view 
we know that any conservation law should reflect a symmetry in the Lagrangian and so far we do not know such a 
symmetry for DE. Thus, it is natural to assume the interaction between the two dark components of the universe. 
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In order to study the interacting model we consider the energy balance equations as 

p m +3Hp m = Q, (22) 
p D +3Hp D {l + w D ) = -Q, (23) 

where Q represents the interaction term which allows the energy exchange between the two dark components of the 
universe. The form of Q is a matter of choice and we take it as 

Q = 3b 2 H{p m + p D ) = 3b 2 Hp D (l + u), (24) 

with b 2 being a coupling constant. Taking the above equations into account and following a same steps as the previous 
section we can find 

WD = -2-n D -^{ 1 + ^-^)- (25) 

When b = 0, wd reduces to its respective relation in the absence of interaction. The first interesting point about the 
equation of state parameter of the GGDE is that in the interacting case independent of the interaction parameter, b, 
for < £ < 1, wd can cross the phantom line in the future while £Id — > 1. 

The deceleration parameter in the presence of an interaction term can be obtained by substituting (|25[) in (|13|) and 
using ([11]). We find [11] 

1 , 3 Q D ( 2b 2 e \ (0R , 

Once again it is clear that setting 6 = 0, respective relation in the previous section is retrieved. Finally, we would like 
to obtain the evolution equation of DE in the presence of interaction. It is a matter of calculation to show that [39| 



a ina 



1-Q D ( 2b 2 £, \ b 2 



2 — — £ V J ^ 



(27) 



Now we consider the stability of this model by study the sign of squared sound speed v 2 . From Eq. (I25|) we have 

dw ■ n 2 D + (2b 2 -0(2Q D -2 + ^ 
W D = 17T7 il D = -RTT^ o- I ^2 Ud - y 2% > 



Also from Eqs. (gSJ) and (JHJ) one finds 

P 1 



P 3H(l + w D + ^Y 



(29) 



Taking into account above relations, 4r = H , and Eq. (f27|) and replacing these relations into (TT8|) it is a matter 
of calculation to show that 

Once again setting £ = leads the results of the original GDE in |28| . The evolution of v 2 against f2 d has plotted 
in Fig. ([2]) for b = 0.15 and £ = 0.1. This figure reveals that v 2 is always negative except for a small period at the 
beginning of the universe and thus, as the previous case, a background filled with the interacting GGDE seems to be 
unstable against perturbation. This implies that we cannot obtain a stable GGDE dominated universe. However, the 
equations show a sensitivity to the value of £ and we found that it is possible to obtain a stable GGDE dominated 
universe but it may be in contrast with theoretical values of /3 [38[ . One important point is also the sensitivity of the 
instability to the coupling parameter b. The larger b, leads to more instability against perturbations. 



IV. STABILITY OF THE INTERACTING GGDE IN NON-FLAT UNIVERSE 

Although it is a general belief that the geometry of the universe is flat, recently this assumption is challenge d by 
several observational evidences. For example, a closed universe is marginally favored by observations from CMB [48j |. 
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FIG. 2: This figure shows evolution of squared sound speed v 2 s versus flo for interacting GGDE model. The solid curve 
corresponds to the GGDE for b — 0.15 and £ = 0.1. Dashed curve is plotted for the original GDE. 



Besides, the measurements of the cubic correction to the luminosity-distance of supernova measurements support the 
idea of a closed universe [49j . Thus, there exist enough motivations to consider DE problem in non-flat universe. 
Here, we would like to extend the study of the GGDE model to a universe with special curvature. In such a case the 
first Friedmann equations read 



H z + — = 



1 



3M| 



(Pm +Pd), 



(31) 



where k is the curvature parameter with k = —1,0, 1 corresponding to open, flat, and closed universes, respectively. 
Taking the energy density parameters ([2]) into account and defining the energy density parameter for the curvature 
term as Qk = k/(a 2 H 2 ), the Friedmann equation can be rewritten in the form 



1 + Slfc = il m + tip. 
Using the above equation the energy density ratio becomes 

Pm fi m 1 + f^fe — SId 



PD 



D 



n, 



(32) 



(33) 



The second Friedmann equation read 



H= -4irG(P + p) + — . 



(34) 



Following the previous sections one can easily see that the equation of state parameter in the interacting non-flat 
GGDE may be obtained as 



1 



w D = 



2 — Or 



2-11 + ^111 + -^ 



2b 2 

o-(i + n*) 



(35) 



From the second Friedmann equation, (|34|) . One can easily obtain 



H 3 

-— r = -fife + -fir>[l + u + w D ] 
H z 2 

and therefore the deceleration parameter in this case is obtained as 



(36) 



H 3 
q = -I - — = -1 - Q k + -Q D [l + u + w D ] 



Substituting Eqs. (|3"3"|) and (|35[) in (|37p we obtain 



q =2 {1 + nk) - 2(2-n D -o 



a 7 



9. k 



2b 2 



SI, 



(37) 



(38) 
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FIG. 3: The evolution of the squared sound speed v\ versus £Id for interacting GGDE model in non-flat FRW universe. The 
solid curve corresponds to the GGDE with b = 0.f5 and £ = O.f , fife = O.Of. Dashed curve shows the original GDE model for 
6 = 0.0.15 . 



Finally the evolution equation of the GGDE in a non-flat interacting case can be written as 



din a 



. - « ■ , , ■ , -(1 + Qk)\ (1 + £lk 

3 2-n D -t\ v *W V 3 y n D K k> ) n D v * 



n k l-firj / A £ \ / fl k \ 26 2 . n A 6 2 



(39) 



In the limiting case fifc = 0, the results of this section,restore their respective equations in flat FRW universe derived 
in the previous sections and for £ = the formulas of [27| are retrieved. The interested reader can see [39| for detailed 
discussion about this case. Here we consider main task of this section, study the stability of interacting GGDE model 
in a universe with spacial curvature. From the energy balance equation (|2"3")l we can obtain 



P = X 

P 3H(l + w D + g-)' 



(40) 



Taking the time derivative of Eq. (l35j) . yields 



n|(i _ £k) + (26 2 (i + n k ) - gi + %)) (2ti D - 2 + p rt 

(2 <>/> • ii-<->7 



^ = - " gv ^ 3/ ' v X;' ^ (4i) 



Using the above relation, replacing ([33)) in Eq. (14T))) , and inserting them into (fT5]l . after using the time derivative 
version of (|39)l one gets 



2 = (n c -i)(e-n D ) 2 (3n c - 4 + 0(i + n fe ) _ 2 (n g + Q 2 - rt g - 3g 
Ws (^-2 + o 2 (n D + o + (o fl _2 + o 2 (fi D + 3 fc (n D -2 + o 2 (n2j + 0" 

Setting f2fc = = b the above relation reduces to the flat noninteracting respective relation. Also the squared sound 
speed of the flat interacting case can be retrieved when Qj- = 0. In order to obtain an insight on the stability issue 
of the interacting GGDE in a non-flat FRW universe, we should consider the sign of v 2 s during the evolution of the 
universe. To this end we plot w 2 versus fin, where the value of SId indicates different epoches of the universe evolution. 
Fig. [3] shows the result which clearly indicates an almost same behavior as the flat interacting case and it seems that 
the presence of the curvature term does not lead to a significant difference in the stability issue at least from this 
feature (the squared sound speed). The squared sound speed, (w 2 ), is positive for < £ and is negative otherwise, 
indicating the instability of the universe against perturbations in the GGDE background. One can easily find from 
(1421) that increasing b will result more instability in the universe. Finally we can say that the stability of a universe 
filled with GGDE and matter against small perturbations crucially depends on the value of £ parameter which is 
defined on the base of the (3 parameter in Eq. ©. 



V. SUMMARY AND DISCUSSION 



Due to the lack of observational evidences every DE model which can explain the current acceleration of the universe 
seems to be acceptable. However, every new model accompanies by many consequences which should be explored. In 
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this paper we considered an important feature, which every DE models should provide. Based on cosmic observations 
our universe is in a stable DE dominated epoch. In this paper we tried to see if the GGDE model is capable to result 
in a stable dominated universe. The main point made the authors to consider such an issue in the GGDE comes 
from the result obtained in [38, 3j|. In these papers authors found that the second term in the GGDE energy density 
([6]) , has a negative contribution to the energy density. This negativity is in contrast of the first term and we thought 
may such a behavior could result in some signs of stability in the model. To this end, we used the squared sound 
speed (Vg = as the main factor for studying the stability. If v 2 s is positive the GGDE would be stable against 

perturbations. When v 2 is negative we encounter the instability in the background spacetime. We have discussed 
several cases including whether there is or not an interaction between dark matter and GGDE and whether there is 
or not a curvature term in the background metric. We found that the stability of the GGDE model crucially depends 
on the parameter £ and adjusting £ to suitable values the model is capable to result a stable DE dominated universe. 
However, observational constraint on the GGDE model [HI rejected mandatory range for £ to leave a chance of stable 
GGDE dominated universe as the present state of the universe. As a result, the universe filled with dark matter and 
GGDE component cannot lead to a stable GGDE dominated universe. We also observed that the instability of the 
interacting GGDE increases with increasing the interacting coupling parameter b. Thus, in comparison to the original 
GDE, the GGDE just has theoretical chance to lead a GGDE dominated universe, however this is ruled out by the 
present observational constraints. 
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